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■ We consider how the measurement of top polarization at the Tevatron can be used 

CN , to characterise and discriminate among different new physics models that have been 

suggested to explain the anomalous top forward-backward asymmetry reported at the 
• Tevatron. This has the advantage of catching the essence of the parity violating effect 

characteristic to the different suggested new physics models. Other observables con- 
C^l structed from these asymmetries are shown to be useful in discriminating between the 

(— I ' models, even after taking into account the statistical errors. Finally, we discuss some 

signals at the 7 TeV LHC. 
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1 Introduction 



The study of the third generation quarks in the Standard Model(SM) fermions continues to 
throw up surprises. Be it a ~ 3o" deviation in A^p^ in the SM fit to the electroweak(EW) pre- 
cision measurements at LEP [1] , a forward-backward asymmetry in top-pair production [2H6] 
significantly larger than what the SM predicts or the discrepancy in the value of the B^-B^ 
mixing as suggested by the recent measurement of the inclusive dimuon asymmetry [7] at the 
Tevatron, some deviations from the SM remain. While, individually, none of them are large 
enough to merit the status of a discovery of physics beyond the SM, they are, nonetheless, 
intriguingly poised to warrant serious attention. This is especially so on account of the large- 
ness of the third generation fermion masses. With the top quark mass scale being very close 
to the EW scale, it is quite conceivable that the third generation (or, at least, the top quark) 
plays an important role in electroweak symmetry breaking itself. Indeed, very many ideas 
going beyond the SM quite often also predict discernible new physics effects for processes 
involving the third generation. On the experimental front, constraints on the universality of 
interactions are far less restrictive for these fermions as opposed to those for the first two 
generations. Hence, on very generic grounds, top quark physics studies at the Tevatron and 
the LHC, hold potential for probing new physics [Sj. It is, therefore, a small wonder that 
very many "explanations" of the physics responsible for the aforementioned discrepancies 
have been offered and a majority of them accord a special role to the third generation. In 
this article, we concentrate on the anomaly in the forward-backward (FB) asymmetry in 
top-pair production. 

Within the SM, the dominant ti production mode at the Tevatron is a pure QCD one and is 
FB symmetric at the tree-level. Weak interactions are highly subdominant and give a very 
small asymmetry {Aps 1%). It is the interference of the tree- level QCD amplitudes and 
higher order terms that yields the dominant SM contribution to this asymmetry and results 
in ApB ^ 5% [Hllin]- This has led to several authors [II1II2] advocating the measurement of 
ApB as a tool for probing physics beyond the SM. 

The CDF and the DO experiments at the Tevatron kindled a great deal of interest in top-pair 
production by reporting a FB asymmetry (Afb) of 17% [2] and 12% [3] respectiveljU. Many 
possible New Physics (NP) scenarios have been offered as explanations. In a later update, 
CDF revised this number to 19.3%. [1] Although the analysis of more data has reduced the 
significance of the CDF result to about 1.8a {Aps = 15%) [3] and DO has presented a newer 
value of 8% [Hj , the deviation from the SM is still non-negligible and, coupled with the other 
longstanding discrepancies in certain third-generation observables, continues to elicit much 
interest [IMH]. 

Forward-backward asymmetry, at the Tevatron, is defined as 

_ ff(cosgt > 0) - ff(cosgf < 0) 

~ aicosOt > 0) + aicosOt < 0) ^ ' ' 

where 6t is the angle made by the top quark with the direction of the proton in the lab-frame. 
Experimentally, the measurement is made in the semi-leptonic channel, where the angle [6^) 
made by the hadronically decaying top (anti-top) with the proton beam and the charge Qi of 

^Thc measurements reported by DO are uncorrected for kinematic acceptance. 
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the decay lepton from the anti-top (top) are together used to construct the net top current 
in the direction of the proton beam, cos 9t above is thus equivalent to —Qi ■ cos 9^ [1] . CP 
invariance is assumed to hold good. 



Naively, a non-zero ApB seems to be an indication of some violation of a discrete symmetry 
and, indeed, most models that purport to explain this anomaly have invoked a parity- 
violating interaction for the top-quark. While this assumption certainly holds true for any 
s-channel NP contribution to qq — )■ ti, clearly, it is not applicable when t- or w-channel 
contributions are present as well. In other words, the measured Aps-, in the presence of any 
NP interactions, may accrue from either explicit parity violation(dynamics) or the effects of 
t-(M)-channel propagators(kinematics) or a combination of both. We believe this issue has 
not been stressed sufficiently in the literature. The new physics explanations have differing 
levels of parity violation encoded in the chiral structure of the interactions. As this is reflected 
in the ApB varying amounts, it would be interesting to construct a probe thereof. The 
polarization of a single top offers one such probe. In fact, this observation was flrst made 
in Ref. [T^ in the context of sfermion exchange contributions to ti production in i?-parity 
violating MSSM (an analog of this model is one of the candidates for an explantaion of the 
observed Apb)- 



It is a thankful coincidence that the top quark polarization is also a quantity that is amenable 
to an experimental measurement due to the large mass of the t. Being heavy, the top 
quark decays before it hadronizes and thus the decay products carry the memory of its 
spin direction. This correlation between the kinematic distribution of the decay products 
and the top spin direction, can be used to get information about the latter. In fact, many 
studies exploring the use of the polarization of the top quark as a probe and discriminator 
of new physics [20] , as a means to sharpen up the signal of new physics [21] and to obtain 
information on tt production mechanism [T9 | l22 y 23] exist in literature. Different probes of the 
top polarisation which use the above mentioned correlation have been constructed [25H2Fj . 
the angular distribution of the decay leptons providing a particularly robust probe due to 
its insensitivity to higher order corrections [2Sj and to possible new physics in the thW 
vertex [27H29]. 



A single-top polarization asymmetry(y4p) can be deflned as 

g(+)-ff(-) , . 

- .(+) + .(-) ■ 

where + or — denote the helicity of the top quark and the helicities of the t are summed 
over. The SM prediction for this arises due to EW effects and is expected to be small. 
In this note, we explore the predictions for Ap originating from the different new physics 
explanations of the FB asymmetry. Spin polarization studies are a part of the agenda at 
the Tevatron as well as the LHC [30]. Both the CDF [31] and the DO [32] experiments at 
the Tevatron have reported measurements of spin correlation coefficients. In addition, CDF 
also reports a measurement of tt helicity fractions [31]. We would like to point out that, 
while these observables may also be able to distinguish between the different NP scenarios 
under consideration, they involve measurement of the polarization of the top as well as the 
anti-top and hence, are experimentally more challenging. This is underscored by the fact 
that the aforementioned measurements are accompanied by fairly large error bars. On the 
other hand, measurement of requires knowledge of the polarization of only the top and 
thus provides an advantage in terms of the statistics that may be obtained. 
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The rest of the article is organised as follows. In the next section, we discuss the rudiments of 
the models we use as templates, with particular emphasis on the features that are germane to 
the issue at hand. This is followed, in section El by a comparative discussion of the structure 
of some observables and their efficacy in distinguishing certain features. The numerical 
results, pertaining to the resolving powers of various observables, are presented in section |H 
Finally, we summarise in section [61 



2 Model Templates 



Rather than probe each and every model that has been proposed, we select some that, to 
our mind, serve as templates. Broadly speaking, four classes of explanations have been 
suggested. The first two involve new t-channel exchanges in qq ti while the third involves 
s-channel contributions to the same. For very large masses of the exchanged bosons, all three 
reduce to four-fermion contact interactions which constitute the fourth category. We do not 
discuss the last-mentioned explicitly as it can be approached from each of the other cases in 
the appropriate limillfl. The particle exchanged in the t-channel could either be a scalar or 
a vector and we shall discuss an example of each. Scalar exchanges in the s-channel do not 
contribute to ApB- On the other hand, s-channel vector exchanges are fairly commonplace 
in new physics scenarios and this constitutes our third template. Finally, we omit tensor 
exchanges because no such well- motivated model leading to large Aps exists. We now discuss 
each template in turn. 



2.1 A flavour-nondiagonal Z' 



While additional U{1) gauge symmetries are well-studied and very often well-motivated, 
it is easy to see that most common variations would not lead to a substantial ApB- In 
models where the Z' appears in the s-channel in the qq — >■ ti process, the NP amplitude 
cannot have a non-zero interference with the QCD contributions. Consequently, obtaining 
the required Aps requires large Z' couplings resulting in too large a correction to cr(tf). On 
the other hand, a fiavour- changing Z' [13] would appear in the t-channel in, say, uu — )■ ti. 
The corresponding contribution io ApB has two sources, viz. kinematic (due to the t-channel 
propagator) and, possibly a dynamic one as well (if the Z' coupling is chiral). Analogous icZ' 
or cuZ' couplings are disfavoured as these would set up fiavour- changing neutral currents. 
Similarly, consistency with i?-physics observables is simpler if the Z' does not couple to the 
6-quark necessitating that the coupling to the top be right-chiral, thus leaving us with a 
Lagrangian of the form [13] 

C3gxZ'^u^^PRt + h.c. . (2.1) 

The model also includes a small fiavor- diagonal coupling to u quarks in order to avoid 
experimental constraints from like-sign top production data from the Tevatron P3] . However, 

■^In principle, there could be scenarios wherein more than one such NP effect could play a role, albeit to 
different degrees. Once again we desist from discussing these explicitly as the gross features thereof can be 
deduced by judiciously combining the templates that we do examine hereafter. 
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this as well as issues regarding anomaly cancellation are not relevant to the discussion at 
hand and, hence, we shall ignore them. 



2.2 Diquarks 



Particles carrying a baryon number of ±2/3 occur in many models, most frequently in those 
concerned with grand unification |33]. While both scalars and vectors are possible, the latter 
would, typically, have a mass as large as the symmetry breaking scale. The scalars can be 
light, though, and may couple to a pair of quarks. Indeed, an example of such couplings 
can be found within the minimal supersymmetric standard model (MSSM) if i?-parity were 
not conserved. Once again, the presence of such couplings would generate a non-zero ApB, 
for both kinematic and, if the coupling is chiral, dynamic reasons. The fermion assignments 
within the SM ensures that the couplings are indeed chiral in nature. Of the various different 
diquark fields that are possible, clearly the one that couples a u to a t can generate the 
maximal amount oi Aps for given mass and coupling strength. In other words, the relevant 
piece of the Lagrangian is [2] 

£ 3 f {ys + yp 75) u + h.c, (2.2) 

where Ta denotes the appropriate color-coupling structure. As can be expected, the cross- 
sections for transforming as a 6 of SU (3)c differs from that for a 3 only in colour factors. 
We will not discuss the two situations separately and will restrict to the color triplet diquarkB. 
While Ref. [H] admits generic ys and yp, it is easy to see that consistency in the 6- sector 
motivates the choice of right-chiral couplings, viz, ys = yp. 



2.3 Axigluons 



Originally motivated as residues of unifiable chiral color models [31], wherein the high en- 
ergy strong interaction gauge group SU{3)l <^ SU{3)p is spontaneously broken to the usual 
SU{3)c, the axigluon was nothing but the octet gauge boson of the broken symmetry 
having a purely axial vector coupling — of strength Qs — with the SM fermions. This model 
has been probed, at Tevatron, both in the dijet |35] as well as the ti channel [T2l[36ll3Z] and 
masses ^ 1 TeV have been ruled out. Although such a scenario immediately predicts App, 
it has the right sign for this quantity only for small values of rriA, well below the Tevatron 
limits [12]. For larger m^, the sign of App flips. 

Motivated by this, Ref. [13] considered a variation with a different embedding of the color 
group into a SU{3)a ® SU{3)b. The gluon and axigluon are now admixtures of the SU{3)a 
and SU{3)b with a non-trivial mixing angle 6a- Anomaly cancellation requires a fourth 
generation of quarks. The new axigluon is 'flavor non-universal' as its couplings with the 3''°' 
and 4*^* generation quarks are different from those with quarks of the first two generations, 
the essential trick being to reverse the sign of the axial coupling of the top-quark, thereby 

''Ref. [M] also investigated scalars (with vanishing baryon numbers) in 8 and 1 representations of SU{3)c 
but those were not found to be very successful in consistently reproducing the relevant experimental data. 
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reversing the sign oi ApB- In the bargain, additional vector- hke couphngs are introduced for 
the axigluon, leading to 

^^^7M(^7y + ^7l75)T,V^A^ (2.3) 

While the vector coupling is generation universal {gy = —Qs cot 26 a), the axial coupling is 
not, with g\ = —Qs cosec26'yi for the first two generations and = +gs cosec26'A for the 
last two. Note that, unlike in cases of the diquark and the Z', wherein the new physics 
contribution only appears in the sub-process uu — )■ ti, for the axigluon, all quark flavors are 
involved. 



3 Analytical Issues 



As we are interested in the polarization asymmetry, we must calculate the cross-sections for 
different final state helicity combinations. While the details are presented in the Appendix, 
note that, in each case, the square of the amplitude can be decomposed as 

\M{Xt, Xt)\' = A+\t\tB + {Xt - Xt) C . 

where Xt, Af = ±1 are twice the helicities of the top (antitop). For observables relating to a 
single polarization, such as that we are interested in, the other must be summed over. 



It is easy to see that while the total tt cross-section and Aps receive contributions from 
only from A, the polarization asymmetry Ap C/A. It is, thus, worthwhile to examine the 
dependence of these terms on the coupling constants. As far as the diquark is concerned 
(see eq. IA.4|) . the terms A (and B) are even in both the (pseudo)-scalar couplings ys (vp), 
and thus, the ensuing App is insensitive to their {ys,p) sign (and, hence, to the chirality 
structure of the theory). The term C being proportional to the product ysVp, the polarization 
asymmetry picks up the chirality structure unerringly. And while eq JA.5[ does not reflect 
this explicitly (owing to the fact that a.V + A structure was chosen), the story is similar for 
the Z'. 



For the axigluon, the situation is bit more complicated. There exist pieces in A (and B), 
that are odd in the individual axigluon couplings, but being parity odd, most of them do 
not contribute to att- The only such piece that contributes to au does so only away from 
the resonance — with the sign reversing as one crosses it — resulting in a subdominant contri- 
bution. Of course, App is very sensitive to the relative signs of the couplings in question, 
which was the original motivation of this model vis-a-vis the flavour-universal axigluon. Once 
again, the sign of would provide additional information about the aforementioned rela- 
tive signs, thereby serving to establish the actual structure of the embedding of SU{3)c in 
the larger gauge group. 



Finally, note that, for each of the three cases, the piece C has a different angular dependence, 
including a parity-odd piece. This leads one to hope that the construction of a rapidity- 
dependent Ap would lend additional resolving power. 
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4 Numerical Results 



From the discussions in the preceding section, one expects, in general, that each of the three 
models would be associated with a different correlation between the three observables a, 
ApB and Ap. However, since the particular values of these quantities depend sensitively on 
the parameters of the model (i.e. the boson masses and couplings), it is conceivable that the 
models could lead to similar values for these observables for some unrelated points in the 
parameter space. Since such a degeneracy would lead to the models being indistinguishable 
from each other, at least as far as such observables are concerned, we begin by delineating 
the situations where this might occur. 

To this end, we scan the parameter space for each of these models with the restriction that 
the couplings be perturbativfl For the Z', this implies (^x < Svr. For the diquark theory, (to 
ensure that there are no relative normalization factors between the couplings in the different 
models under consideration) the same limit is imposed on gy = \/2y where, y = {y's + v'pY^'^ ■ 
In the case of axigluons, the perturbativity condition is imposed on the couplings qa and qb 
(associated with SU{?>)a and SU{3)b, respectively) with both qa, Qb < Stt. This, in turn, 
translates to 10°< 9 a < 45° [I5]. 

As far as masses are concerned, the lower limits are model-specific. It is easy to appreciate 
though that, for very large masses, all the models would, essentially, be equivalent to contact 
interactions. For very large masses, engendering a large enough ApB would, then, need very 
large couplings, thereby coming into conflict with the perturbativity requirement discussed 
above. In effect, only M„e^ ^ 3000 GeV is relevant. Furthermore, we impose direct limits 
by taking into account various considerations as detailed below. 

In case of diquarks, one could argue for m$ > 350 GeV, motivated by the limits on squark 
searches at the Tevatron [T]. However, for a scalar decaying primarily into a top (or a 
bottom) and a jet, these limits do not apply and the corresponding limits are expected 
to be much weaker. Although such a dedicated study is yet to be reported, preliminary 
investigations |38] suggest that m$ ^200 GeV cannot yet be ruled out, and we shall adopt 
this. 

As for the Z' model, much depends on the flavor structure of its coupling. If the Z' is 
very light, then the branching fraction for t — )■ Z' + n is no longer negligible. Further, if 
the Z' decays hadronically, this would lead to an enhancement in the contribution to the tt 
production from the lepton+jets channel as compared to the contribution from the dilepton 
channel giving rise to deviations from the experimental measurements in the two channels. 
Ref. [13] has already considered this to find that rriz' ^120 GeV is consistent with current 
data and we shall adopt this limit as well. 

The limits on the flavor non-universal axigluon mass need more thought though. Exper- 
imental searches, so far, have only considered flavor universal axigluons. Based on these, 

^We impose this condition only until the 10 TeV scale. Were we to demand that the theory remain 
perturbative until a significantly higher scale, the parameter space would be further restricted. However, in 
the absence of any knowledge of the ultraviolet completion, we desist from this. 
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ttla ^ 1250 GeV is ruled out by direct searches while the distortion in the mu spectrum [39] 
is too severe for ^ 1400 GeV. However, it should be noted that the resonant production 
cross sections for flavor non-universal axigluons differ on account of the differences in their 
coupling to light quarks. Nonetheless, we do restrict to niA ^ 1400 GeV. 



The parameter spaces thus defined, the numerical calculation is done using a parton level 
Monte Carlo routine. We use the CTEQ6L [IQ] parton distributions with the factorization 
scale set to rrit with the later held to 172.5 GeV to be consistent with the value used in 
measurements of cross-section [H] and Aps [H]- A i^T- factor [12] of 1.33 is used to estimate 
the cross- section at NLO|^. 



< 




-0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 -0.05 0.05 0.1 0.15 0.2 0.25 0.3 0.35 

(a) (b) 

Figure 1: (a) Correlation between Ap and App for different models. The vertical solid (dot- 
ted) lines correspond to the central value and 1-a bands of the new (old) CDF measurement 
of ApB, namely IS-ffVo ± 5.5% (19.3% ± 6.9%). The 'star' corresponds to the SM value at 
NLO. For each model, the scan is limited to only perturbative coupling strengths, (b) As in 
the previous panel, but restricted to only those points that are consistent with the experimen- 
tally observed cross-section at the 1-a level and with restrictions on Mboson as described in 
the text. 



For each parameter-space point we calculate au, App and as defined above. Our calcu- 
lation of App includes the SM part [S]. As for Ap, the SM contribution is miniscule. This 
is primarily because a non-zero Ap can arise only when the production mechanism treats 
positive and negative helicity states differently. Consequently, a pure QCD process can never 
give rise to a non-zero Ap. In the SM, Ap will arise only due to electroweak effects, either 
at the tree level or through mixed EW-QCD higher order corrections, and is, hence, rather 
small (C(10~'^)). Given this, we restrict ourselves to the tree-level value of Ap as indicated 
in the previous section. 



In Fig 1(a), we show the correlation between Ap and App for each of the three models 
alongwith the experimentally allowed 1-a band for App. No restrictions have been imposed 
on the physical observables, barring the aforementioned perturbative limit on the couplings. 
Note that the splayed-finger like structure for $ is but a reflection of the sampling density 
in the masses and a finer sampling would have resulted in a dense region in this plane. It 



^In the absence of NLO calculations that incorporate the new physics effects under consideration, we use 
the SM K-Factor. 



8 



immediately transpires that the three models under discussion have very different correla- 
tions between the observables. However, this is evidently misleading as a large faction of the 
parameter space depicted here gives rise to a. ApB well outside the experimentally allowed 
range. This is even more so for the total cross section. Clearly, large values of the asymme- 
tries would occur only for substantially large NP couplings and these, typically, give rise to 
too large a cross section. 



Imposing the restriction that the total ti production cross section should lie within la of the 
experimentally observed value and restricting M^oson as described earlier, eliminates bulk of 
the parameter points and leads to the correlation plot of Fig, 1(b) One feature immediately 
stands out. For each of the models, the allowed regions split into two, one close to the SM 
point and the other distinctly separate. 
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Figure 2: Variation in a(tt) and Afb with coupling for different choices of Mi;,. 

For the Z' and the diquark, the reason for this is easy to understand. Clearly, for small values 
of the NP coupling, the cross-section is close to that within the SM and, hence, consistent 
with the measured value. However, for such small couplings, the ApB generated is also small 
and not in consonance with the experimental observation. For a given boson mass, as one 
considers larger values of the coupling, sufficient Afb can be produced while the cross-section 
is prevented from becoming too large because of the destructive interference between the SM 
and NP pieces. In fact, for an intermediate range of the coupling, the cross section falls well 
below the SM value. The location and the extent of the range naturally depends on the 
mass of the boson. For the case of the diquark, this is illustrated in Figj2l The situation is 
analogous for the Z'. 



For the axigluon, the situation is somewhat more complicated. Here the interference may be 
constructive or destructive depending upon the values of s and Ma under consideration. For 
Ma > 1400 GeV, the interference is mostly destructive though. For 9a ~ 45°, the situation 
is identical to that in the case of the flavor universal axigluon, with the exception that now 
the relative minus sign between the g\ and g^^ yields ApB that is of the right sign as well 
as the right magnitude. As 6a decreases, the destructive interference takes over leading to 
a{tt) values that are too low (see Fig, 3(a)). For Ma < 2000 GeV, a further increase in the 
strength of the couplings allows the pure NP term to become dominant and the cross-section 
increases once again to become consistent with data. Due to this, a part of the range of 6a 
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below 25° gets ruled out. Again, the actual location and extent of the range depend on the 
particular value of M^. On the other hand, for 6*^^ 20°, ApB is only produced in small 
amounts (Fig, 3(b) ). The resultant picture that emerges is the tail of points near the bottom 



left of Fig 1(b) which are consistent with the measurement of a{tt) but not with that of ApB- 
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Figure 3: Variation in cross-section and ApB with 9 for different choices of Ma- 



It is interesting to note that, in Fig, 1(b) , for the regions that are consistent with the measured 
ApB values (whether the old measurement or the new one), there is no overlap between the 
three models. It is true that, consistently flipping the relative sign between the axial- 
and vector-like coupling^, would flip the sign of while a and ApB remain unaltered. 
However, even after such a variation, the flavor non-universal axigluon case would still retain 
its distance, in the ApB ~ Ap plane, from the other two models. Thus, a measurement of Ap 
should, in principle, be able to distinguish the models even before the corresponding boson 
has been directly observed, this being particularly true of the axigluon. 



Let us also note at this point that analyses of Ref. [23] indicate that polarization values of 
around 15%-20% can be measured at the Tevatron with 10 fb~^, with good significance. 
Further, in addition to the azimuthal asymmetry explored in Ref. [23], a simple forward 
backward asymmetry of the decay lepton w.r.t. the beam direction, may also be used at the 



Tevatron and hence one would be able to increase the sensitivity further. From Fig. 1(b) 
one sees thus that the values of polarization predicted by the different models seem large 
enough to afford measurement at the Tevatron. 



While discrimination of the axigluon model from the other two seems clear, that between 
the diquark and the Z' models is not so straightforward. For one, they lead to overlapping 
ranges in both Ap and App- Furthermore, FigsH] have been drawn with the error bars 
suppressed. Incorporating the latter would render the two apparently disparate population 
zones overlapping. If, however, one were to flip the relative sign between the scalar and 
pseudoscalar couplings of the diquark, the dissimilarity between the models would become 
quite conspicuous. 

^Of course, this would necessitate a different embedding in the gauge group. 
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Figure 4: Ra vs A for forward-backward asymmetry and polarization asymmetry. The points 
depicted are consistent with the CDF measurementc of a{ti) and ApB o,t the 1-a level. Also 
imposed are the perturbativity constraint on the couplings and the constraints on Mboson 
described in the text. 



Nevertheless, even without changing the relative signs of ys and yp and with a view to making 
a more robust distinction between the two models, we propose observables that exploit the 
rapidity distributions in the two scenarios. With the typical values of the t-channel masses 
being different in the two scenarios (see Figs. 2 & 3), it is but natural that the typical rapidities 
would be different. To this end, we define the ratios Rap and Rapb where 

_ A{\Ay\ < 1) 
^ A{\Ay\ > 1) ' 

and Ay is the difference between the rapidities of the top and the anti-top. The CDF 
Collaboration does report a measurement of Aps in different Ay regions [5], leading to 
R^^^ = 0.043±0.194 with the SM expectation being Rapb{SM) = 0.317±0.053. While this 
seems to offer a potentially strong discriminator, as we shall soon see, the associated errors 
do not, yet, allow for discrimination between models. Instead, we advocate an examination 
of correlations involving Rap and Rapb- 

In FigJH Rapb Rap are plotted against App and Ap respectively. For consistent com- 
parison with CDF results [S], the asymmetries are calculated in the ti rest frame and include 
the NLO contribution from the SM [S]. It can be seen that the diquark and Z' models 
populate different regions in the Rapb ^ ^fb and Rap — Ap planes and the separation is 
quite distinct. 

Once again. Fig. |4] does not incorporate errors in the measurement of the quantities plotted. 
The statistical part thereof can be easily estimated. We obtain the efficiencies (i.e. includ- 
ing kinematic cuts, event selection efficiencies etc.) from the number of events reported 
in the CDF App measurement [S], and assume that the corresponding efficiencies for Ap 
measurement would be similar. Calculating the expected number of events in this way, we 
then add the errors in quadrature. While this might seem a trifle optimistic, we compensate 
by restricting our projections to an integrated luminosity (~ 4.2 fb~^) less than what has 
already been analysed in Ref. [5]. 
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(a) (b) 

Figure 5: Correlations involving Afb, Ra^b' ^'^^ -^^p o,ft^i" inclusion of error bars. The 
points depicted are consistent with the same constraints as explained in the caption of Fig!^ 



In FigEl we redraw FigJH along with the error bars computed as described above. Although 
the central values in the two scenarios are widely separated, a significant overlap between 
the error bars persists, thus makes it difficult to distinguish between them. Indeed, the size 
of the errorbars also shows that the Rapb measurement, on its own, is, as yet, incapable of 
conclusively distinguishing the models from even the SM. On the other hand, the correlation 
between and Rafbi (FigE]) looks more promising. While some overlap exists still, the 
extent of overlap is distinctly smaller. While a quantitative analysis confirms this, we desist 
from sharpening our conclusions at this stage. It should also be realized that an analysis 
more sophisticated than what we have performed and including all the data on tape, is quite 
likely to further reduce the relative sizes of the errors. 



1.1 
1 
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^ 0.6 
< 

^ 0.5 
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0.1 


0.1 0.2 0.3 0.4 0.5 0.6 0.7 

"Afb 

Figure 6: Correlation between Rafb ^'^^ Rap- The points depicted are consistent with the 
same constraints as explained in the caption of Fig^ 

It should be realized that Rap is not the only direction-dependent asymmetry variable. Of 
the many such possible, we consider only one namely as calculated separately for the 
forward and backward hemispheres. In particular, we plot, in FiglTJ the ratio of in the 
two hemispheres against the total Ap. Again, one finds that the diquark and Z' models give 
rise to markedly different correlations. 
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Figure 7: Correlation between the ratio of Ap in the forward and backward hemispheres and 
total Ap. The points depicted are consistent with the same constraints as explained in the 
caption of Fig\^ 



However, as the experience with R^p has shown, the effect of the separation is hkely to be 
diluted by the errors involved. With experimental errors under control, a combination of 
observables as described above, can be expected to distinguish between the "explanations" 
for the ApB "anomaly" quite successfully. 



5 Observables at the LHC 



The initial state being symmetric at the LHC implies that no simple forward-backward 
asymmetry w.r.t the beam direction can be defined. Hence, we limit ourselves to the cor- 
relation between Ap and (T(ttjl. For the parameter space points that are consistent with 
the Tevatron measurements of cross-section and App a.t the 1-a level, we find that the three 
models correspond to different correlations (see Fig. [8]). In going from the Tevatron to the 
LHC, the relative contribution from the gg — >■ tt sub-process (which does not contribute to 
such asymmetries) increases as compared to the qq — )■ ti subprocess. The quark initiated 
subprocess suffers a further suppression as anti-quark fluxes are diminished at a pp collider. 
As a result, the magnitudes of at the LHC are lower than those observable at the Teva- 
tron. Nonetheless, it is interesting to see that the separation of the expected values of 
into different islands is a trend that persists. That the Ap for the axigluon is now almost 
universally negative (as oposed to the Tevatron where positive Ap were allowed) is but a 
consequence of the fact that a resonance is more easily hit at the LHC. 

It should also be borne in mind here that the values of polarization presented here correspond 
to those obtained by integrating over the entire kinematical range. It is clear that by making 
suitable cuts on the observables such as rriu and/or pip, the polarization may be enhanced 
substantially, as has been seen in Ref . [23] . It is also clear that the optimal cuts will be again 
different for different production mechanisms of the NP top pairs. 

^The K-Factor in this case is estimated by taking the ratio of the NLO cross-section quoted in Ref. [43] 
and our own LO calculation resulting in K=1.8. 
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Figure 8: Correlation between Ap and a{ti) at the LHC {^/s = 7 TeV). The vertical lines 
show the 1-a interval of the CMS J^3| / (dot-dashed) and ATLAS l44^ (dashed) measurements 
of the ti cross-section. The points depicted are also consistent with the constraints mentioned 
in the caption of Fig^ 

At the LHC, there also exists the possibihty of direct detection of these models through 
the production of the corresponding boson. For diquarks and axigluons (albeit, the flavor 
universal variety), the production cross-sections have been calculated in Refs. [H] and [12] 
respectively. For the Z\ direct detection is possible through pair production of Z's or 
production of a Z' in association with a t. We present the corresponding cross-sections 
for these processes in Fig. [9l assuming that eq. fl2.1[) is the only term in the interaction 
Lagrangian. 




[TeV] Mz, [TeV] [TeV] 

Figure 9: Several production cross sections at the LHC for the Z' model. Note that obtaining 
the requisite ApB at the Tevatron requires the coupling to be typically larger (gx ^ Ij- 

Although the cross-sections for the first two processes seem quite large, claims about de- 
tectability should be made cautiously. Note that such Z' would decay into a top and a light 
quark (unless other couplings are switched on). Thus, these production processes would, 
essentially, lead to a tt pair accompanied by one or two hard jets. The QCD background for 
the same is quite large. While invariant mass reconstruction would, in principle, distinguish 
between signal and background, in practice this is not an easy task , especially with two 
top-quarks present. This observations holds for di quark production as well. 
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Relatively easier to look for is the production of a like-sign top-pair, for which this model 
predicts a large rate. Indeed, as Fig. M^c) shows, with the projected luminosities at the LHC, 
such a Z' could be measured upto several TeVs even for gx = 1 (which, essentially, takes 
us to the limit of a contact interaction). It should be noted, though, that a very large mz' 
necessitates a large gx for producing the correct ApB- Such a scenario should have made its 
presence felt in the very same process at the Tevatron itself. 



6 Summary 



The anomalously large forward-backward asmymmetry in top-pair production observed at 
the Tevatron continues to puzzle. Several models have been proposed to "explain" this. 
While some of them have their roots in well-motivated scenarios originally proposed to 
address other issues, some of the others are purely phenomenological in nature. The very 
fact that many different models of new physics can explain this anomaly renders difficult the 
identification of the best solution. 

In this note we have analysed the role that the longitudinal top polarization (Ap) can play 
in discriminating between such models. The scenarios proposed, typically, differ from each 
other in their chiral structure and in the relative amount of FB asymmetry that is generated 
by the dynamics as opposed to kinematic effects. Ap, being a pure parity violating effect, 
probes the former aspect. 

To illustrate this point we choose from among the different proposed models, three templates 
and calculate the predictions for in each. As expected, we find different correlations 
between Ap and ApB in the three cases and some scenarios can be clearly distinguished 
from others purely on this account. However, bearing in mind that experimental errors 
may reduce this apparent separation between the models, we further construct the ratios 
Ra = A{\Ay\ < l)/A{\Ay\ > 1), with Ay being the rapidity separation between the top 
pair. Correlation between — Rafb seems quite promising in its ability to separate the 
different models even after the inclusion of errors. Given that such measurements are already 
being made at the Tevatron, we hope that the adoption of our algorithm would serve to solve 
this vexing issue. 

We also look at the case of the LHC briefly. Again, a scan over the parameter ranges of 
the three template models, consistent with the Tevatron data, shows that they occupy three 
disjoint islands in the Ap — ati space. Although the magnitude of is smaller, owing to 
the increased importance of the gg — )■ ti subprocess, the situation could be improved further 
by the adoption of apropriate kinematic restrictions. In this particular exploratory exercise, 
as far as the LHC is concerned, we refrain from carrying out such studies. 

While most of the models to explain App had been chosen by the respective authors, so 
as to evade current direct search bounds, this would not be the state of the matter in 
the LHC era. As we have demonstrated, the production cross-sections for new particles 
inherent to such models, are likely to be large at the LHC. Given the peculiar couplings 
of such particles, commensurate search strategies need to be devised, but, by no means. 
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is this an insurmountable problem. Indeed, even if the new particles were all heavy, the 
effective Lagrangian responsible for generating the Aps would, typically, leave its mark on 
SM processes at the LHC (see SecJS]). 

It should be noted that, by their very nature, most of the scenarios proposed are phenomenog- 
ical in nature and have concentrated on explaining Aps while taking care of consistency with 
only the most obvious of other observables. Indeed, a global fit would render many of them in- 
consistent. For example, the diquark model would take Re = T[Z — > hadrons)/r(Z — )■ 
away from the LEP/SLC measurements |15]. As for the Z' model, the couplings and masses 
required to give a good fit to both au and the new measurements of are disfavoured by 
the lack of observation of like-sign top pairs at the Tevatron [13]. The flavor non-universal ax- 
igluon, on the other hand, is disfavoured by 5-physics observables [TTj. While such problems 
might seem specific to these models, variations of these and others such are almost endemic 
to all such endeavours. They can, of course, be cured by introducing compensatory effects in 
the respective models. Thus, rather than rejecting them on account of such disagreements, 
direct searches are the best bet. 

It is here that algorithms such as those we propose are particularly useful. Apart from 
discriminating betweeen models, such analyses can also serve to form an information basis 
on which more sophisticated model-building can be based. 

Note Added : As we were finalising the manuscript, two papers, [IH] and [17], have 
appeared. In the former, the authors have studied different polarization observables, one of 
them being the polarization of the t quark that we consider. While they look at different 
models proposed to explain the Aps-, they concentrate solely on the LHC (at 7 TeV as well 
as 14 TeV) and do not consider correlation with ApB- Ref. [17], on the other hand, looks at 
the effects at the Tevatron but it is an extension of model-independent analysis of Aps done 
in their earlier paper [18] to predictions for correlations of t and i polarizations. Since CP 
conservation relates t and i polarizations, we consider it sufficient to study the polarization 
of either one of t or t, which gives us the advantage of larger statistics. Further, the analysis 
of Ref. [im is valid only when masses of exchanged particles are much higher than the cut 
off scale and hence there is no direct, easy comparison between the two approaches. 

After the submission of our manuscript, the CDF Collaboration published a new analysis of 
ApB [IS] including its dependence on the mass and rapidity of the ti pair. The analysis in 
different rapidity regions is identical to that in Ref. [5] which we have already considered. 
Further, Ref. [H] reported that ApB found to be larger in the high invariant mass region 
('^tf ^ 450 GeV) than in the low invariant mass region (m^j < 450 GeV). While Ref. [H] 
itself pointed out that the analysis is still fraught with theoretical uncertainties, we carried 
out a preliminary investigation which showed that this feature is borne out by all the model 
templates that we have studied (Fig JTOj) and does not offer any additional resolution between 
the models under consideration. In fact, this was to be expected in view of the facts that (a) 
the total cross section, which is dominated by the small mti range agrees very well with the 
SM expectations and (&) the measured ApB is substantial. Furthermore, even if uncertainties, 
both theoretical and experimental, can be reduced, the crucial fact is that, unlike Ap, this 
observable has little sensitivity to the chiral nature of the new physics couplings. 
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Figure 10: ApB in the regions rriu < 450 GeV and mu > 450 GeV compared to the total 
ApB- ApB is calculated in the tt rest frame. The points depicted are consistent with the 
same constraints as explained in the caption of Fig\^ 
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Appendix 



We present here the squared matrix elements for the process qq — >■ ti, where, in the sub- 
process center of mass frame, t has a velocity /3 and subtends an angle 9 with q. The final 
state quarks carry polarizations \t and Af respectively. These can then be used to define the 
requisite density matrices. For each case, one can write 



\M\^=\MsM? 



NP 



[k.l] 



where the first part on the r.h.s. represents the SM contribution while the second represents 
the additional contribution due to the new physics (including the interference with the SM 
amplitude, and hence not positive definite). For the SM piece, we have. 



\M 



SM\ 



91 



(1 - \t\t){l + cl) + (1 + \t\i) (1 - (5^) (1 - cl) 



(A.2) 



where cq = cos 9 and sq = sin 9. 



In expressing the lA^l^rp for the respective cases, it is useful to introduce the notation 



T = l-(3ce 



I (ml - 1) 



U = l + (3ce 



(A.3) 



Triplet Diquark Exchange 



With y"^ = Vs + Up, "we have 



\M\%p 



192 (u - MlY 



y'^U'-A yi. y'p Xt Xt (/? + c,)^ + 2 ys yp y' (A, - Af) f/ (/3 + c,) 



+ 



9s s 



72 (u - Ml) 



y'i:^-P'ci) 



-Xt Xiy^ (2 cl + si) + 2ysyp {Xt - Xt) (/3 + 2cb + ^c]) 
For the sextet diquark, only the color factors would change [T 



(A.4) 
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Z' Exchange 



Defining x = / 2 m|, , we liave 



\M\%p 



|f/2 + 2x(l-/32)+x2r2} 
XtXi {{l3 + ce)^ + 2x{l-l3^)cl + x^^-ce)^} 



16 {t - M|,)2 

+ {\t - \t) + ce)U + 2 X (1 - c, - x^ {(3 - ce)T} 



9s9x s 



18 {t - Ml) 



{U^ + l-/3' + x (r2 + l-/32)} 

\t\t{il+x) (2 4 + (3^ sl) + 2il-x)(3ce}_ 
+i\t - Af) {/3 (1 - x) (1 + cl) + 2 (1 + x) ce} 



Flavor Non-Universal Axigluon Exchange 

Defining combinations of couplings as 

we have 



\M\ 



Q 1 

9v 9 a 
9v 9a 



sJa 



'Pv = 9v9v 
Pa = 9Wa 



NP 



18 [{s - MIY + T\M\] 



where 



Ja 



SqS+ (1 + P^cD+SqS^ (1 - + ^VQVTPce 
-Xth [Sq S+ {/3' + cl)+SQSr{l- 13') c] + SVqVt (3 ce] 



+ 9l{s-M\) 



{AVA(3ce + 2rv{2-(3^ si)} 

- At Af {4 Pce + 2Vv{2cl + f s^) } 
+ (At - \i)WAg\,ce + 2glg\f3{\ + c^) 



For brevity's sake, we are not presenting the SM expressions for gg — )■ tt. 
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